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We investigate the effects of the nearly fulfilled Efimov conditions on the properties of three-body 
resonances. Using the hyper-spheric adiabatic expansion method we compute energy distributions 
of fragments in a three-body decay of a nuclear resonance. As a realistic example we investigate the 
1~ state in the halo nucleus 11 Li within a three-body 9 Li + n + n model. Characteristic features 
appear as sharp peaks in the energy distributions. Their origin, as in the Efimov effect, is in the 
large two-body s-wave scattering lengths between the pairs of fragments. 

PACS numbers: 21.45.+V, 31.15.Ja, 25.70.Ef 



Introduction. The Efimov effect was introduced more 
than thirty years ago as an anomaly in a three-body sys- 
tem arising when at least two of the three two-body s- 
wave scattering lengths approach infinity [l| . Then an in- 
creasing number of three-body bound states appear close 
to the two-body threshold even if there are no two-body 
bound states. The effect is prohibited by the Coulomb 
potential while only diminished by higher angular mo- 
mentum Q. Although entirely possible in molecules Q 
the effect is unlikely to appear in nuclei due to unfavor- 
able mass ratio 0, |5| . 

Still there exists a number of nuclear systems, called 
halos 0, which are natural three-body systems - a core 
plus two neutrons - where the Efimov condition of at 
least two large scattering lengths is nearly fulfilled. Al- 
though the unfortunate combination of the heavy core 
and light neutrons prohibits the appearance of bound 
Efimov states in the discrete spectrum, they still may 
appear as peculiar structures in the continuum. Very 
little, however, is known theoretically about the Efimov 
states in the continuum. 

Experimentally, on the other hand, the number of ac- 
curate and kinematically complete experiments for three- 
body decays of nuclear resonances is rapidly increasing 
0, • Also three-body decays of excited states of small 
molecules are presently experimentally investigated in 
details . The measured observables are the width and, 
particularly promising, the energy distributions of the 
three fragments after the decay Although a num- 
ber of theoretical studies involve calculations of three- 
body continuum properties [ToL HTl IT^| , calculations of 
energy distributions for resonances under Efimov condi- 
tions have not been done before. In contrast to bound 
states, investigations of the fingerprints of the Efimov ef- 
fect on the decays of three-body resonances are so far 
lacking. 

In this letter we report on an investigation of nuclear 
three-body resonances under the nearly fulfilled Efimov 
conditions. In particular we calculate the energy distri- 
butions of the decay fragments of a nuclear resonance 
and trace the origin of the characteristic peaks in these 



distributions. 

Calculation of energy distributions. We assume that 
three particles emerge after decay of a preformed res- 
onance state. At large distance we then strictly have a 
three-body problem. This is not necessarily true at small 
distance where the three-body cluster assumption may 
be inappropriate. We shall extend the concept from two- 
body nuclear a-decay. There the outer part of the po- 
tential between the daughter nucleus and the a-particle 
is known and the inner part is adjusted to give the cor- 
rect resonance energy. This treatment accounts for the 
major variations of a-decay widths. The fine-tuning is 
obtained by using the preformation factor describing the 
probability for finding an a-particle at the inner turning 
point of the two-body potential. For two-body decays 
the width is determined by the outer part of the wave- 
function, while the fragment energy is fixed from energy 
conservation. It is anticipated and intuitively plausible 
that for three-body decays both the width and the frag- 
ment energy distributions are determined by the large 
distance behavior of the wave- function. 

The notion of large distance is not a priori obvious for 
three particles where either all three or only two inter- 
particle distances can be large. We shall specify distances 
by the value of the hyper-radius p, 
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where m; is the mass and the cm. coordinate of the 
particle number i, and m is an arbitrary mass scale. The 
other five hyper-spheric coordinates are dimensionless an- 
gles, f2, which determine the directions and relative val- 
ues of the coordinates of the constituents . 

Within the hyper-spheric adiabatic method one distin- 
guishes the fast, J7, and slow, p, coordinates. Then, for 
every fixed slow coordinate p the eigen-value problem is 
solved for the fast angular coordinates O 



(2) 



where H(p) is is the Hamiltonian of the three body sys- 
tem with fixed p, W n {p) are the angular eigen- values, and 
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Q) is the angular eigen- function with a parametric 
dependence on p. In practice instead of the Schrodinger 
equation J5J) we solve the mathematically equivalent Fad- 
deev equations |l3| . 

Then $ n (p, fi) are used as a basis for the total wave- 
function ^>(p, f2), 



(3) 



where the expansion coefficients satisfy the system 

of hyper-radial equations where the eigen- values W n (p) 
of the fast angular subsystem serve as effective potentials, 



£l + ^ + ^(W n (p)-E) ]./„(,! 
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where £7 is the three-body energy and P nn ' are the non- 
adiabatic terms @|- 

A three-body resonance corresponds to a (complex- 
energy) solution ^(p, 0) of the system (0} with the 
asymptotic boundary condition of an outgoing wave in 
every channel n, 



f n {p — > oo) = C„ exp(+irep) 
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where C„ is an asymptotic normalisation coefficient and 
k = yj2mE/h 2 is the three-body momentum, which is 
the conjugate of p. The other five variables in momentum 
space are the momentum angles ft K which are conjugates 
of the coordinate space angle fl. The angles ft K determine 
the directions and relative values of the momenta of the 
fragments. In practice instead of the exponential function 
in we use a corresponding Hankel function. We also 
employ the so called complex scaling |13( with an angle 
6, p — > e p, which improves the numerical accuracy of 
the calculations. 

An example of the effective potentials and the hyper- 
radial resonance function is shown on Fig. ^ The lowest 
potential has an attractive pocket where a resonance re- 
sides. The corresponding dominating component of the 
hyper-radial function shows a typical resonance behavior 
- a pronounced maximum in the pocket region with sub- 
sequent oscillations, matching the corresponding Hankel 
function, in the asymptotic region. 

As shown in [l4| the momentum distributions P(fl K ) of 
the fragments can be calculated from the wave-function 
ft) as 
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For the graphical presentation we shall integrate out the 
four directional angles in fl K and only leave one angle 
which specifies the energy Ei of the fragment i relative 
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FIG. 1: The lowest adiabatic potentials W„(p) for the 
11 Li(l~) halo nucleus within the three-body 9 Li+n+n model 
with interactions from [l5| . The n-core and n-n scattering 
lengths are a nc ~ a nn ~20 fm. The inset shows the low- 
est hyper-radial resonance function with its large distance 
asymptotics - the Hankel function. The complex scaling an- 
gle 9 — 0.15, the resonance angle 9r — 0.12 corresponds to 
the resonance energy of about 0.4 — O.li MeV. 



to the maximum possible energy ^ max ' of this fragment 
in the decay process. 

Different decay mechanisms carry distinct signatures 
of the process, e.g. sequential decay via a two-body reso- 
nance must give a peak in the energy distribution of the 
third particle at the peak-energy preferred by the first 
two-body decay. However, when the two-body s-wave 
interaction is only weakly attractive - no bound states 
exist - no stable intermediate structures are available for 
such a sequential decay. Yet some characteristic effects 
might appear especially when the threshold for binding 
is approached for more than one two-body s-state, i.e. 
when the Efimov conditions are nearly fulfilled. 

Investigated systems. The nuclear candidates for the 
Efimov conditions are systems of two neutrons and an 
ordinary core-nucleus where as needed only short-range 
interactions are present. These structures are found as 
ground or excited two-neutron halo states at the neutron 
drip-line or for more stable nuclei 0- 

We take the (1~) state in the 11 Li halo nucleus within 
the three-body 9 Li+n+n model as a realistic example. 
The interactions from |15| provide the n-core and n-n 
scattering lengths of about 20 fm and the effective ranges 
of about 5 fm. 

The adiabatic potentials W n (p) for this system are 
shown on Fig.^ The lowest potential shows an attractive 
pocket with a relatively long tail, extending up to about 
20 fm, that is p < a nci despite the fact that the binary 
interactions are short-range Gaussians with the range of 
about 5 fm. This long tail is due to interference between 
the subsystems with long s-wave scattering lengths 0. 
This tail would lead to the bound Efimov states if only 
the mass ratio of the constituents were inverse, or if the 
scattering lengths were larger. 

We solve the hyper-radial equations with the 



3 



boundary condition JSJ and obtain the resonance hyper- 
radial wave-function. The dominating lowest component 
is shown in the inset of Fig. ^ 

From the hyper-radial wave-function we extract the 
asymptotic normalization coefficients in and calcu- 
late the energy distributions of the core, 9 Li, and the 
neutrons using the expression 

In the realistic calculation apart from the large scat- 
tering length s-wave interactions also the interactions 
in higher partial waves, including spin-orbit and tensor 
forces are included. Those interactions might affect the 
distributions and mask the fingerprints of the Efimov 
conditions. We therefore in addition present two model 
calculations - one with only s-wave neutron-core inter- 
actions, and one where also the s-wave neutron-neutron 
interaction is included. Thus only effects of the large 
scattering lengths can be seen in the distributions. 

For the model systems we restore the resonance energy 
by increasing the scattering length by a factor of two 
(a nc pa a nn pa 50 fin) and by tuning the three-body force. 
The latter is short-range and therefore does not affect the 
asymptotics of the system. The larger scattering lengths 
also serve the purpose of amplifying the visible signatures 
of the Efimov conditions. 

Convergence. The limit p — > oo in (j^J is reached in 
the region where the tail of the lowest adiabatic potential 
becomes negligible, that is somewhere about and beyond 
a nc . In practice we calculate the distributions at some 
finite Pmax ~ o-nc and then check that an increase in p max 
does not alter the distributions. 

In the hyper-spheric adiabatic method the increase in 
Pmax demands a corresponding increase in the number 
N of the adiabatic potentials which have to be included 
in the wave-function. Thus the convergence has to be 
checked for the two parameters, p max and N. 

In our numerical convergence tests we calculated the 
distributions with parameters 3 < N < 12 and a nc J$ 
Pmax ~ 2a„ c and checked that there is a region in the 
parameter space, where the distributions are unchanged. 

The model systems converge much easier and demand 
smaller p max and N. 

Numerical results. Let us first look at the model sys- 
tem with only s-wave n-core interactions. The energy 
distributions of fragments are shown on Fig. [21 where E 
is the fragment's energy and I?( max ) is the largest possi- 
ble energy. The convergence for this system is reached 
fairly well already at p m ax = «nc and N = 2. 

The core distribution shows one sharp peak at 
E c /Ec ~ 0.5 while the neutron distribution shows 
two sharp peaks, at E n /E n ma ^ pa and E n /En pa 1. 
Since only s-wave interactions are allowed these peaks 
must be the sought fingerprints of the large scattering 
lengths or, in other words, virtual states. 

The geometric interpretation of these peaks is that the 
decay proceeds by emission of one neutron with maxi- 
mum energy, 



£(max) = E 
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FIG. 2: The energy distributions of the fragments - the core, 
9 Li, and the neutrons - in the decay of a three-body res- 
onance 11 Li(l _ ) calculated in the three-body 9 Li + n + n 
model with only s-wave n-core interactions (scattering length 
a nc pa 50 fm). The different curves are calculated with dif- 
ferent pmax and different numbers of adiabatic channels N to 
illustrate the convergence. 



(where E r is the energy of the resonance, m c and m„ are 
the core and the neutron masses) allowing the exploita- 
tion of the other neutron-core virtual state with vanishing 
relative energy. This produces the neutron peak close to 



E n l E n 



(max) 



1. 



The other neutron and the core together pick the rest 



of the energy, E r — E n ma " 



E r 



and share it 



proportional to their masses, since they have vanishing 
relative energy and thus the same velocity. Consequently 

and the other 



the core gets E, 
neutron gets E n 



-- E r 

E r 



m c -\-m n -\-m n m c +m ri 



Since m n <C m c we have for the predominant core 



energy 



E c — E r 



m. 



where 



m c + m n + m n m c + m r , 
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Thus the peak in the core distribution should be, in 
agreement with calculation, just below E c /Ec = 0.5. 
The other neutron has the predominant energy 



E n — E r 



m c + m n + m„ m n + m c 



< £( max > (10) 



which results, again in agreement with the calculation, 
in a second peak around E n /E n ' 



(max) 



0. 



The neutron-neutron s-wave interaction will allow an 
additional configuration, the two neutrons with vanish- 
ing relative energy against the core with maximum en- 
ergy. Consequently, another peak should appear in the 
£ c /£^ max) distribution close to E c /E^ ] pa 1. 

The two neutrons get the rest of the energy, E — 



E, 



(max) 



E, 



r rn c -\-rn n -\-m r 



which they share equally. Thus 



4 



3 

X) 



— i — i — i — 1 — i— 

Li-core 

p m u -1.0 a : N=3 

p = 1.5 a m ; N=5 

p =2.0 a : N=5 



— i — 1 — i — 1 — r~ 

neutron 

p = 1.0 a ;N=3 

p = 1.5 a ;N=5 

p =2.0 a ; N=5 



J_ 



J_ 



0,25 0,5 



0,75 



E /E 



,(max) 



0,25 0,5 0,75 1 

E /£ (max) 

n n 



FIG. 3: The same as fig. |3|with additional s-wave interaction 
in the n — n subsystem (a nn » a nc ~ 50 fm) . 
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FIG. 4: The energy distributions of the fragments - the core, 
9 Li, and the neutrons - in the decay of a three-body resonance 
11 Li(l~) calculated with interactions from fl5| where the n- 
core and n-n scattering lengths are about 20fm. 



a neutron will predominantly have the energy 
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-E r 
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m c + m n + m n 
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which again should result in an additional peak just be- 
low £„M max) « 5- 



These additional structures are clearly seen in Fig. |3| 
where the corresponding numerical distributions are 
shown. Indeed there are two additional peaks at 
E c /Ec w 1 and E n /En ~ |- Convergence is 
somewhat harder to achieve compared to the first model 
system and one can see that for p max = 2a nc five adia- 
batic channels are not quite enough. 

If the virtual states are the dominating properties of 
the system, the full realistic calculation should exhibit 
structures, similar to Fig. [3] And so it does, as seen on 
Fig. 0] where the realistic distributions are shown. De- 
spite the peaks being somewhat blurred due to interfer- 
ence with higher partial waves, the characteristic struc- 



tures at E c /E[ 
and En/E^ 11 * 
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E c I E c 



(max) 
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I can still be distinguished in the real- 
istic distributions and can thus serve as an experimental 
indication of the nearly fulfilled Efimov conditions. 

The convergence for the realistic system is the hard- 
est to achieve numerically - one has to go up to p max = 
2.5a nn and consequently include N = 12 adiabatic chan- 
nels to get the converged distributions. This calculation 
includes about 960 different hyper-spheric partial waves 
in the basis for the wave-function and seems to be con- 
verged within about 5% accuracy. 

Conclusions. We have investigated the effects of the 
nearly fulfilled Efimov conditions on the properties of 
three-body resonances. We have calculated the energy 
distributions of fragments in a three-body decay of a res- 
onance for a realistic system, the (1~) state in the halo 
nucleus 11 Li within a three-body 9 Li + n + n model, and 
two model systems where only s-wave interactions were 
allowed. We have shown that characteristic features ap- 
pear as sharp peaks in the energy distributions and have 
traced their origin to the large scattering lengths between 
the decay fragments, like in the Efimov effect. The ex- 
perimental verification of the predicted distributions is 
straightforward but perhaps challenging due to the re- 
quired accuracy. 
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